Conductance oscillations in zigzag platinum chains 
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Using first principles simulations we perform a detailed study of the structural, electronic and transport prop- 
erties of monoatomic platinum chains, sandwiched between platinum electrodes. First, we demonstrate that the 
most stable atomic configuration corresponds to a zigzag arrangement that gradually straightens as the chains 
are stretched. Secondly, we find that the averaged conductance shows slight parity oscillations with the number 
n of atoms in the chain. Additionally, the conductance of chains of fixed n oscillates as the end atoms are pulled 
apart, due to the gradual closing and opening of conductance channels as the chain straightens. 

PACS numbers: 73.63.-b,73.40.Jn,68.65.-k,71.15.Ap 



The existence of single atom chains was demonstrated 
some time ago using the scanning tunneling micro- 
scope (STM) and mechanically-controllable break junctions 
(MCBJ) JlJ where a quantized conductance close to Go = 
2e 2 /h for gold was measured, in agreement with previous 
theoretical predictions |2j|. Since then, a number of exper- 
iments 0] and theoretical calculations ^ have proved 
that the 5d elements Au, Ir and Pt can be used to produce 
monoatomic chains. For gold chains 1 1], the average distance 
between conductance peaks was found to be 2.5 A with the 
conductance G of the last plateau being very close to Go- The 
ensemble-averaged conductance also shows small oscillations 
around Go as the length of the chain increases |4]. 

For platinum, the average distance between the peaks in 
the length-histograms is about 1.9-2.3 A In contrast 

with gold, the conductance is no longer an integer multiple 
of Go, but instead decreases from 1.6 Go to 1.2 Go as the 
length of the chain increases jH. In addition, significant con- 
ductance oscillations are superimposed on top of this decreas- 
ing trend, which have been attributed to a parity effect. 0,0. 
The rich behavior of platinum compared with gold arises from 
the larger number of conduction channels at the Fermi energy 
{Ef) due to the presence of d bands|7]. 

To understand the structural, electronic and transport prop- 
erties of platinum chains, we have performed a complete se- 
ries of first principles simulations of platinum chains attached 
to platinum fee leads. We have employed our newly developed 
code SMEAGOLlH %, which calculates the density matrix 
and the transmission coefficients of a two probe device us- 
ing the non-equilibrium Green's Function formalism (NEGF) 
H 1 Oil . The scattering potential is calculated self-consistently 
using the SIESTA implementation of density functional the- 
ory Illlfl2lfl3ll . We have approximated the exchange and cor- 
relation potential by the Local Density Approximation (LDA), 
since we have found that provides a slightly more accurate de- 
scription of the structural and conducting properties of bulk 
and infinite platinum chains than GGA 11411 . We have not 
included in our calculations spin polarization which is only 
relevant for very stretched chains 111 511 - Our main result is that 
the most stable arrangement of platinum chains corresponds 



1 1 1 1 
(a) 


\ ^ 

v ^_ ~ — 


- (b) 

1 


' 1 ' 1 ' 1 ' 1 ' - 


1 , 1 , 1 , r 



I 1 I 1 I 1 I 1 I 1 L_ 

,4 1,6 1,8 2 p 2,2 2,4 
d (A) 




E-E p (eV) 



FIG. 1: (color online) (a) Cohesive energy E c and (b) number of 
open scattering channels at the Fermi energy N(Ef) of platinum 
chains as a function of d z . Number of open scattering channels 
N(E) as a function of energy E calculated at (c) d z = 2.38 and (d) 
d z — 2.15 A for zigzag (solid line) and linear (dashed line) chains. 



to zigzag configurations, that are straightened as the chains 
are stretched. The conductance decreases from 1.6 Go to 1.2 
Go upon chain stretching, showing slight parity oscillations. 
Aditionally, we also find geometry-induced oscillations as the 
electrodes are pulled apart. Such oscillations are due to the 
gradual opening and subsequent closing of transport channels 
as chains with a fixed number of atoms n are straightened. 

We start by discussing the case of infinite platinum chains, 
which is useful to understand the basic physics of the con- 
striction when a long chain is formed. We use a two atom unit 
cell with periodic boundary conditions along the three spatial 
directions. The cell is large enough along the x and y axes 
to avoid spurious interactions between the chain in the central 
unit cell and its images. We simulate two kinds of chains: lin- 
ear chains, where the atoms are constrained to lie along the 
z axis, and zigzag chains, where forces are allowed to relax 
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along the three spatial coordinates. We find that zigzag chains 
are more stable than linear chains, as it can be seen from the 
cohesive energy of Fig. ^ a )- The equilibrium distances along 
the z axis are d z eq = 2.15 and 2.38 A for the zigzag and lin- 
ear chains, respectively. In the zigzag arrangement the Pt-Pt 
bonds are located in the xz plane and make a 24.8 degrees an- 
gle with the z axis. This means that the interatomic distance 
for the zigzag configuration is d = 2.37 A, which is almost the 
same as for the linear case. These angles decrease almost lin- 
early as the chains are stretched, and become approximately 
zero for d z larger than 2.5 A. Interestingly, if the distance is 
reduced below the zigzag minimum the system falls into an- 
other stable configuration with a ladder arrangement, similar 
to that predicted by Sen et al. 1 16]. 

Fig- Hb) shows the number of open scattering channels 
N(Ep) at the Fermi energy Ep of an infinite zigzag chain, as 
a function of d z . The figure reveals that N(Ep) = 5 for small 
d z , but decreases to N(Ep) = 2 as the chain is stretched be- 
yond 1.95 A, and maintains that value for quite a large range 
of distances. For d z somewhat larger than 2.30 A, N(Ep) in- 
creases in two steps to 4 Go, and stays constant thereafter until 
the chain breaks. To understand which channels are opening 
and closing, we have studied the projected density of states 
(PDOS) for chains at different stretching. As limiting cases 
we compared a linear chain at d z = 2.38 with a zigzag chain 
at d z =2.15 A. For the linear chain we find that the d xy and 
d x 2_ y 2 orbitals are completely filled, while the hybridized s, 
d X z, d yz and d z 2_ r 2 orbitals have all finite weight at Ep, lead- 
ing to four open channels, as shown in Fig. ^c). Atoms in 
zigzag chains at d z = 2.38 A make small angles with the z- 
axis at this distance and therefore, the d z -dependence of the 
PDOS and N(Ep) look both fairly similar to those of the lin- 
ear chain. In contrast, the number of open channels N(E) as 
a function of energy E, shown in Fig. ^d), and the PDOS of 
zigzag chains at d z = 2.15 A look rather different from those 
of the linear chain. The zigzag chain has only two, and not 
four, open channels at the Fermi energy, corresponding to a 
mixture of all d-orbitals, since now the s-orbital is completely 
filled, while the d^j, and d x 2_ y 2 have moved up in energy. 
This analysis of infinite chains can explain why in some ex- 
periments the conductance increases when the electrodes sep- 
arate, because there is a gradual transition from a zigzag to 
a linear configuration upon the stretching of the chain, with 
the linear chain presenting larger number of open scattering 
channels. 

This simplified picture cannot explain however the features 
that are present in a real experiment, where the contact to 
the electrodes and rearrangement of the leads near the sur- 
face can decrease the transmission of channels or even close 
them. To investigate this possibility, we have simulated finite 
chains of various lengths (between 1 and 5 atoms), attached 
to fee platinum leads oriented along the (001) direction (see 
Fig. 13. The leads are composed of repeated slices of 3 x 3 
atoms and are connected to the chain through a square of 4 
atoms, i.e. through the fee (001) hollow site. In order to get 
rid of undesirable oscillations in the transmission coefficients, 
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FIG. 2: (color online) The different chains connecting (001) oriented 
fee leads studied in this paper: (a) Four-atoms chain, (b) single atom 
contact and (c) 2-atom chain. d z is the distance between leads. 



we have used periodic boundary conditions along the xy plane 
and summed over 12 fc -points. To treat the contact region self- 
consistently, three atomic planes of bulk platinum are included 
in each side of the scattering region| 17]. To calculate the most 
stable configuration for each chain, we perform the relaxation 
by keeping fixed the bulk Pt leads and relaxing the apexes of 
the point contact. In this way the relaxation is performed over 
the chains and the two 2x2 planes forming the hollow site. 
The ground state energy is therefore calculated as a function 
of the distance d z between the outer slices (unrelaxed) as in- 
dicated in figure|2] In Fig.|5]we plot the cohesion curves for 
zigzag atomic chains with a number of atoms ranging between 
2 and 5. 

The ground state energy of the single atom contact (not 
shown) as a function of distance is a parabola whose mini- 
mum is located at the equilibrium distance d z e<i — 6.3 A. A 
parabolic dependence is also found when we place two atoms 
facing each other, with d z eq = 9.1 A. However, if we allow 
the outer planes to move along the x and y directions, these 
two atoms achieve a zigzag configuration that has a lower en- 
ergy, as shown in Fig. [3] This zigzag arrangement is actually a 
local minimum, as in the case of chains of infinite length. For 
short distances, the atoms at the apex gain energy by forming 
further chemical bonds to other atoms at the electrodes. The 
same kind of curves are observed for 3, 4 and 5 atoms, with 
rf 2 ,eq = H-4, 13.3 and 15.4 A, respectively. We find that the 
difference in lengths between chains of n and n + 1 atoms lies 
in the range 1.9 - 2.1 A, in very good agreement with experi- 
ments 1 4]. The only configuration which does not fit into this 
pattern is the single atom contact, whose equilibrium distance 
is 2.8 A away from the two-atom chain. Indeed, the single 
atom contact does not constitute a chain by itself, since a large 
number of bonds link the central atom to its neighbors, which 
must be broken simultaneously in order to snap the chain. In 
contrast, atomic chains can be topologically characterized as 
broken by cutting just one single bond, which leads us to as- 
cribe the peaks found in length histograms |4] to chains 2-, 3-, 
4- and 5-atoms long. Fig. [2] also indicates that the region of 
stability of the (n + 1)— atom chain begins at a distance where 
the the n— atom chain is very stretched and therefore close to 
being broken. This fact helps understand why the first two 
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FIG. 3: (color online) Cohesion curves of the chain plus leads sys- 
tems, for zigzag and linear chains of two atoms (circles and squares, 
respectively), and zigzag chains containing 3, 4 and 5 atoms (trian- 
gles, diamonds and stars, respectively). All curves have been shifted 
in energy in order to align the local minima. 



peaks found in length-histograms, that correspond to chains 
with two or three atoms, are much higher than those attributed 
to chains of four or five atoms 1 3, 4], because there is a signif- 
icant probability that an n-atom chain will break, before the 
n + 1-atom chain forms. At a more detailed level, we also find 
that the angles between the atomic bonds and the chain axis 
are small for short chains (9.1 degrees for the 3-atom chain), 
but increase with the chain length (21.6 degrees for 4-atom 
chain). The case of the 5-atom chain is more complicated: 
the two atoms joining the leads make angles of 36.3 degrees, 
while those in the middle have angles equal to 16.1 degrees. 

Moving to transport, we observe that the single atom con- 
tact behaves very differently from any other, with a conduc- 
tance of about 5.9 Go at the equilibrium distance, which de- 
creases almost linearly to G — 3.5 Go until the contact breaks. 
These values are very different from those of reference j(J, 
where the pyramid-like structure of the contact was not in- 
cluded, probably leading to an overestimate of the interatomic 
distance. For the 2-atom chain we find a conductance of 1.5 
Go at the equilibrium distance in the zigzag arrangement and 
a conductance ranging between 1.8 and 2.0 Go in the linear 
case, as shown in Fig. |4ja). The second value can be as- 
sociated with the return conductance measured when a con- 
tact is made again after the system breaks |6]. We therefore 
propose that the first configuration established between two 
atoms should be linear when both tips are brought together. 

For larger chains, we obtain a non-monotonic behavior of 
the conductance as a function of d z . When the separation be- 
tween the leads increases from a compressed configuration, 
the conductance initially decreases and then exhibits a plateau 
at around the equilibrium distance. Under further expansion it 
grows again (for stretching of about 1 A beyond the equilib- 
rium separation) and finally decreases exponentially. An ex- 
ample is shown in Fig. 0](b), for n = 3. This behavior can be 
understood as follows. For small distances, atoms in the chain 
are very close to each other and there is a large number of open 
channels. As the distance increases, the transmission through 
many of these channels is reduced and the conductance de- 
creases until a value between 1 .0 and 1 .5 Go is reached, where 




FIG. 4: (color online) Conductance of platinum chains, (a) 2 atoms 
in zigzag or linear arrangement (circles or squares, respectively), (b) 
3 atoms; in the inset is plotted the tail of the curve in a logarithmic 
scale, that shows the tunneling behavior when the chain breaks. 
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FIG. 5: (color online) (a) Evolution of the conductance of 2-, 3-, 4- 
and 5-atoms chains as a function of d z (circles, triangles, diamonds 
and stars, respectively), (b) Conductance vs. the number of atoms 
in the chain measured at the equilibrium distance (crosses) and aver- 
aged over a range of two Angstrom about such distance (squares). 

the chain has a zigzag configuration. If the chain is stretched 
further, the conductance increases again, following the evolu- 
tion from a zigzag to a linear chain, until a small plateau with 
a conductance of 2 Go is formed. Finally, the chain enters 
the exponential tunneling regime, as shown in the inset to Fig. 
E](b), whose onset signals the point where the chain would 
break or a new atom would enter the chain. 

To try to make contact with the experiments of Ref. U the 
conductances of all chains are gathered together in Fig. 5 (a). 
We find clear oscillations with a periodicity equal to the inter- 
atomic distance, which have a structural origin. These oscilla- 
tions are due to the gradual closing and opening of channels, 
which occurs as the angles between the atomic bonds and the 
z-axis increase and decrease. For a stretched chain, the angles 
increase if a new atom enters the chain, and subsequently de- 
crease as the chain is further stretched. Moreover the size of 
the conductance is linked to the biggest angle subtended by 
atoms in the middle of the chain. Larger chains tend to have 
bigger angles and lower conductances. One exception is the 
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5 -atom chain, whose conductance, 1.2 Go, is slightly bigger 
than that of the four-atom chain (1.0 Go). This is also easily 
explained by our calculations, since atoms in the middle of the 
5-atoms chains have a smaller angle than those of four atoms 
chains (16.1° versus 21.6°). 

The large geometry oscillations of Fig. 5(a) may mask 
those due to the parity effect, which have a larger periodic- 
ity but a smaller amplitude. The experimental conductance 
curves shown in Ref. 0] are obtained by averaging individ- 
ual conductance traces over thousands of pulling circles. We 
therefore expect that this averaging should smooth out short- 
distance fluctuations, such as the large geometry oscillations 
discussed above, and uncover parity oscillations, that have a 
longer wavelength. To explore such a possibility, we have 
plotted in Fig. 5 (b) the conductance of each chain, taken 
at its equilibrium distance. We find that the large geometry 
oscillations have been completely washed out, leaving only a 
small oscillation that has a periodicity equal to twice the inter- 
atomic distance. To further check the stability of these results, 
we have also plotted in Fig. 5 (b) the conductance of each 
chain, but now averaged over a range of lengths of about two 
Angstroms around each equilibrium distance. We believe that 
such a procedure may reproduce the most important features 
of the experimental averaging of conductance curves. Inter- 
estingly, both curves almost overlap. The remaining curve 
actually shows two peaks, corresponding to chains with three 
and five atoms. We note that the position of the conductance 
peaks, the overall shape of the curve and the magnitude of the 
conductance are in excellent quantitative agreement with the 
results of Ref. U In view of the above results, we propose 
that the geometry oscillations seen in our simulations are av- 
eraged over in the experiments by Smit et al. |4], leaving only 
the smaller-amplitude, but robust parity oscillations. 

In summary, we have shown that monoatomic platinum 
chains have a zigzag structure, both in the case of perfect infi- 
nite chains, where the conductance is halved compared to the 
linear configuration, and the case of chains between two fee 
(001) electrodes. We find that atomic chains are composed of 
at least two atoms, while the single atom contact shows sig- 
nificantly different structural and transport features. Our cal- 
culations explain the negative slope in the conductance that is 
found for larger chains, and accurately reproduce the oscilla- 
tions given by parity effects found in the experiments of Ref. 
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